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ABSTRACT 

Aims. We present mm line and continuum observations at high spatial resolution characterizing the physical and chemical properties of the 
young massive star-forming region IRAS 18182-1433. 

Methods. The region was observed with the Submillimeter Array in the 1.3 mm band. The data are complemented with short-spacing 
information from single-dish CO(2-l) observations. SiO(l-O) data from the VLA are added to the analysis. 

Results. Multiple massive outflows emanate from the mm continuum peak. The CO(2-l) data reveal a quadrupolar outflow system consisting 
of two outflows inclined by ~90°. One outflow exhibits a cone-like red-shifted morphology with a jet-like blue-shifted counterpart where a blue 
counter-cone can only be tentatively identified. The SiO(l-O) data suggest the presence of a third outflow. Analyzing the '^CO/"CO line ratios 
indicates decreasing CO line opacities with increasing velocities. Although we observe a multiple outflow system, the mm continuum peak 
remains single-peaked at the given spatial resolution (~ 13500 AU). The other seven detected molecular species - also high-density tracers like 
CH3CN, CH3OH, HCOOCH3 - are all ~l-2" offset from the mm continuum peak, but spatially associated with a strong molecular outflow peak 
and a cm emission feature indicative of a thermal jet. This spatial displacement between the molecular lines and the mm continuum emission 
could be either due to an unresolved sub-source at the position of the cm feature, or the outflow/jet itself alters the chemistry of the core 
enhancing the molecular abundances toward that region. A temperature estimate based on the CH3CN(12i^ - lU) lines suggests temperatures 
of the order 150 K. A velocity analysis of the high-density tracing molecules reveals that at the given spatial resolution none of them shows 
any coherent velocity structure which would be consistent with a rotating disk. We discuss this lack of rotation signatures and attribute 
it to intrinsic difflculties to observationally isolate massive accretion disks from the surrounding dense gas envelopes and the molecular outflows. 
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1. Introduction 

The formation of massive stars is an important topic by it- 
self, but it is also essential for understanding cluster for- 
mation and the Galactic evolution. Significant observational 
progress has been made over the last decade by statistical 
studies of large samples of young massive star-forming re- 



gions(e.g..iriume e 
2000; "Hunter et al." 
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'2000'; 'Sridharan et al."2002'; 'Beuther et al.' 
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Mueller et al.-,2002- Shirlev et al.-2003 


; .Walsh et al.. 


20031 iFaundez et alJl2004 iHill et alJl2005l iKlein et alJl2005l). 



However, to really understand the physical properties at the 
inner center of the massive star-forming regions at very 
early evolutionary stages, high-spatial-resolution studies at 
(sub)mm wavelengths are a more adequate way of investi- 
gation. Observational studies of various massive star-forming 
regions with different interferometers have started disentan- 



gling the c omplex nature of high-mass star formation (e^ 
Cesaroni et al. 1997: Zhang et al. 1998: Wvrowski et al 199S 



Shepherd et alJl2000l: iBeuther et alj ' l2005bt ICiibb et alJl2003[ 

Garav et al ."2003VS u et alJl2004 . Nevertheless, we are still at 
the beginning of such studies, and there is urgent need to in- 
crease the number of investigated sources at high spatial and 
spectral resolution to base the conclusions on a more statistical 
foundation. To mention a few specific questions we are inter- 
ested to tackle: What is the nature of massive molecular out- 
flows, in an evolutionary sense as well as between sources of 
different mass and luminosity? What is the core mass distribu- 
tion at the earliest evolutionary stages, can we infer information 
about the formation of the Initial Mass Function? What is the 
nature of potentially embedded accretion disks, and which way 
to best investigate them? Here, we present Submillimeter Array 
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(SMA') 1.3 mm continuum and line observations toward the 
young high-mass star-forming region IRAS 18182-1433 ad- 
dressing some of these questions. 

The source IRAS 18182-1433 is part of an intensely 
studied sample of 69 potential High-Mass Protostellar 
Objects (HMPOs) at early evolutionary s tages prior to 
forming a significant U CHii region (Sridh aran etalJ l20Q2t 
iBeuther et alJ l2Q02alblch . At a near kinematic distance of 
4.5 kpc, the luminosity and gas mass from the region are 
~ 10'*^ Lq and ~ ISOO Mp, respectively ( Sridharan et al. 
l20Q2HBeuther et alJl2002al) . The region harbors a bipolar out- 
flow with an outflow mass of ~ 30 Mo and an outflow en- 
ergy of ~ 2.3 X 10'*'' erg ( Beuther et al. 2002b). Furthermore, 
IRAS 18182-14 33 shows OH, Class ii CHc,OH, and H?0 
maser emission ("Forst er & Caswel]lll999t fWalsh et alJll998t 
[Beuther et al . 2002c) but was not detected in cm continuum 
emission down to 1 mJy, indicating th e early evolutionary stage 
of the source jSridharan et alJEo02l) . However, recent sensi- 
tive VLA observations at cm and 7 mm wavelengths detected 
two weak cm and one 7 mm continuum source in this region 
jZapata et alj2005 l). A spectral index analysis of the three sub- 
sources shows that the 7 mm emission likely traces the main 
massive protostellar object in the field, whereas one cm feature 
is consistent with a thermal jet ~2" offset from the 7 mm peak. 
The other cm feature ~ 1 1 " south-east of the 7 mrt i peak shows 
signa tures of non-the rmal synchrotron emiss ion dZaoata et al] 
I2OO5 I. see also Fig.l2ll. lDe Buizer et allhOQ ^ observed this re- 
gion in the mid-infrared at 10.5 //m and 18.1 fim wavelength. 
They detected a weak mid-infrared source close to the 7 mm 
and the cm thermal jet feature. Furthermore, they detect an- 
other far stronger and extended mid-infrared source ~10" to- 
ward the south-east which appears to be associate d with the cm 
contin uum emission with negative spectral index dZapata et alJ 
l2005h . 

2. Observations 

2.1. Submillimeter Array Observations 

The HMPO IRAS 18182-1433 was observed with the 
Submillimeter Array (SMA) on April 30th 2004 with 6 an- 
tennas at 1.3 mm in the compact configuration with pro- 
jected baselines between 9.8 and 53k/l. The phase center was 
R.A.=18''21"'09.098', decl.=-14°31'49".29 (J2000.0) with a 
tuning frequency of 230.55 GHz in the upper sideband (visr = 
59.1 kms '). Parts of the correlator were mis-functioning, nar- 
rowing down the usable bandwidth to approximately 1.6 GHz, 
the covered spectral range was 219.45 to 221.07 GHz in the 
lower sideband (Isb) and 229.18 to 230.72 GHz in the up- 
per sideband (usb), respectively. The spectral resolution was 
0.8 125 kHz/channel corresponding to ~ 1.1 kms ' velocity 
resolution. The weather was excellent with zenith opacities 
t(225GHz) between 0.1 and 0.07 measured by the National 
Radio Astronomy Observatory (NRAO) tipping radiometer 

' The Submillimeter Array is a joint project between the 
Smithsonian Astrophysical Observatory and the Academia Sinica 
Institute of Astronomy and Astrophysics, and is funded by the 
Smithsonian Institution and the Academia Sinica. 



operated by the Caltech Submillimeter Observatory (CSO). 
Passband and flux calibration were derived from Uranus obser- 
vations, and the flux density scale is estimated to be accurate 
within 20%. Phase and amplitude were calibrated with regu- 
larly interleaved observations of the quasar NRAO530 (11.8° 
from the source). We applied diff'erent weightings for the con- 
tinuum and line data resulting in synthesized beams between 
3.6"x2.4" and3.9"x2.5". Therms noise ofthe line and contin- 
uum images is dominated by the side-lobes and the correspond- 
ing inadequate cleaning. The dirty beam with side-lobes at ap- 
proximately a 40% level is shown in Fig.Q] The 3cr rms of the 
1.3 mm continuum image is ~ 15 mJy. The varying rms of the 
diff'erent line images is discussed at the appropriate places in 
sections l3T2l & l3 .31 The initial flagging and calibration was done 
with the IDL superset MIR originally developed for the Owens 
Valley Radio Observatory ( Scoville et al. 1993) and adapted for 
the SMA^ . The imaging an d data analysis were conducted in 
MIRIAD (ISaultetal Jl995l). 

2.2. Short spacings from the IRAM 30 m 

To complement the CO(2-l) observations with the missing 
short spacing information, we combined the SMA data with 
the previously published CO(2-l) map from the IRAM 30m 
telescope ( Beuther et al. 2002b). For observing details of the 
single-d ish observations we refer to the paper bv lBeuther et alJ 
( l2002lj ). The single-dish data were converted to visibilities 
within the MIRIAD package (using the task UVMODEL). 
The single-dish and interferometer data were then subsequently 
processed together. The weighting of the two datasets was cho- 
sen to recover large-scale emission but maintain at the same 
time the high-spatial resolution to resolve small-scale struc- 
ture. Convolving the merged dataset to the 11" beam of the 
30 m data alone, the peak fluxes in the channel map are at ap- 
proximately a 70% level of the single-dish measured fluxes. 
Measuring the integrated fluxes toward the same outflow re- 
gions for the merged dataset and the single-dish data alone, we 
even recover in the merged dataset -88% and -89% of the blue 
and red CO(2-l) single-dish emission, respectively. The result- 
ing dirty beam from the merged dataset is shown in Fig. [2 as 
well. The synthesized beam of the combined data is 4.2" x2.8" 
(PA. 15°). 

3. Results 

3.1. Millimeter continuum emission 

Figure |2presents the synthesized 1.3 mm continuum image of 
the region produced by averaging the apparently line-free parts 
of the upper and lower sideband spectra presented in Figure 
|3l We identify only one massive isolated source at the given 
resolution and sens itivity. The 7 mm continuum source from 
IZapata et al.l ( l2005l) is associated with the 1.3 mm peak posi- 
tion, whereas the cm source with positive spectral index - thus 
likel y a thermal jet featu re - is approximately 2" to the south- 
east ('Zap ata et alJ l200.^ and appears to be associated with 

^ The MIR cookbook by Charlie Qi can be found at 
^http://cfa-www.harvard.edu/~cqi/mircook.html 



Beuther et al.: SMA observations of the HMPO IRAS 18182-1433 



3 





















Q 








VJ 








1 ? .n 



Table 1. Millimeter Continuum Data 
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Fig. 1. The dirty beams of the SMA data alone (top panel) and 
the merged SMA+30m data (bottom panel). The contour levels 
cover the range from +1 in steps of +0.125. 



the elongation of the 1.3 mm emission in this di rection. The 
CH3OH Class 11 maser feature jWalsh et alJll99 8») is probably 
rather associated with the cm jet feature than with the mm con- 
tinuum peak. H2O maser emissi on shows various mas er fea- 
tures distributed over the region jBeuther et alJl20()23) . some 
of them appear to be associated with the cm feature whereas 
others are found closer to the mr n continuum pea k. The weak 
mid-infrared feature detected by IPe Buizer et alJ ll2005 ) may 
be more closely associated with the cm feature than with the 
mm peak. However, the given mid-infrared astrometry (~1.5") 
is not good enough to clearly associate that feature with one or 
the other peak, and we refrain from further interpretation of the 
mid-infrared and mm/cm continuum emission. Although the 
spatial resolution of the 1.3 mm continuum data is insufficient, 
the continuum extension in the direction of the cm jet feature 



S peak 
S int 

Mass@43K 
A'h,@43K 
Mass@150K 
jVH,@150K 



426 mjy/beam 
546 mjy 
47.6 Mo 
5.7 lO^^cm-^ 
12.4 Mo 
1.5 10^3 



indicates that the mm continuum emission is not only of pro- 
tostellar nature but that there can be additional mm continuum 
emi ssion due to an un derlying jet (for other similar examples 
see iGueth et alJl2003l or Beuther etal. 2004b). Nevertheless, 
the majority of the mm continuum emission is due to cold dust 
from the protostellar object(s). 



i[Jy] 

0.4 




n 
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09:0 
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Fig.2. 1.3mm continuum image of IRAS 18182-1433. The 
contours range from 10 to 90% (step 10%) from the peak flux 
shown in Tabled The synthesized beam (3.6" x 2..4"P.A. 
13.9°) is shown at the bottom-left. The star marks the 7 mm 
source and the open and filled triangles the cm sources with 
positive and negative spectral index, respectively dZapata et alJ 
12005). The square and white crosses mark the positions from 
the C H^OH Class n and H^ O masers, respectively (Walsh et alJ 
ll998tlBeuther et all2002cl) . 



Compared with the single-dish flux at 1 .2 mm wavelength 
of 3.9 Jy (^euther et al. 2002a), the integrated interferometer 
flux of ~ 546 mJy is rather low, of the order 13% from the 
single-dish measurements. Since these interferometer observa- 
tions are not sensitive to spatial scales > 25", we filter out 
a rather smoothly distributed large-scale halo, and only the 
compact core remains. Assuming optically thin dust emission, 
we can estimate the mass and column density of this central 
1. 3 mm dust a nd gas co re following the equa tions outlined 
in'Hildeb randljl983.) and iBeuther et al.1 ll2002al) . We calculate 
these parameters assuming two temperatures, 43 K as derived 
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Fig. 3. Upper and lower sideband spectra vector- averaged over 
all baselines. The molecular lines are marked. 



from two-components grey-body fits (IS ridharanet'ai]l2002l) . 
and 150K consistent with the CH3CN(12i - 11^) emission 
(see ^3.51 . Based on previous interferometric observations of 
similar sources from the original sample (IRAS 18089-1732 
& IRAS 20293+3952, Beuther etal. 2004a b), we use a dust 
opacity index y6 of 1 , corresponding to a dust opacity per unit 
dust mass of /c(1.3mm) ~ 1.8cm^g"'. We note that /3 = 
1.4 corresponds to A-(1.3mm ) ~ 0.9cm^g"', as suggested by 
lOssenkopf & HenningI 19941 and using a larger k as we do re- 
sults in lower mass estimates. Depending on the temperature, 
the central core mass is between 12.4 and 47.6 Mq. The H2 
column density varies between 1.5 x 10^^ and 5.7 x lO^"* cm"^ 
(Tabled, corresponding to visual extinctions of a few hundred 
(Ay - A^H2/0.94 X 10^', Frerking et al. 19 82). At a temperature 
of 43 K, the 3cr rms of 15 mJy corresponds to a mass sensitivity 

of - 1.3 Mq. 

We tried to derive radial intensity profiles, but since the core 
is only marginally resolved and the data suffer strongly from 
the missing short spacings, the resulting intensity distribution 
is unreasonably steep prohibiting a further analysis in this di- 
rection. 



3.2. Spectral line emission 

Figure |3] presents the observed spectral bandpass around 220 
and 230 GHz in the Isb and usb, respectively. Altogether we de- 
tect 25 Hnes from 8 species ('^CO, HNCO, H^^CO, SO, ^'*S02, 
CH3OH, HCOOCH3, CH3CN) and 2 additional CO isotopo- 
logues ('^'CO and C"*0) with lower energy levels fiiower be- 
tween 5 and 568 K (Table|2J. The broad range of energy levels 
shows the power of broad spectral bandpasses which allow to 
study very different excitation regimes simultaneously. 

Figure 0] presents integrated line images of all species and 
isotopologues (except ■''^SOa which is blended with CH3OH 
and too weak for imaging). For '^CO and '^CO we are show- 
ing the blue and red line wing emission identifying a bipolar 
outflow centered on the mm continuum peak. While the blue 
'^CO emission is confined to a small region south-east of the 
mm-peak, the red '^CO emission shows cone-like emission to- 
ward the north-west. The outflow will be discussed in detail 
in ^3.31 All other spectral lines show compact emission close 
to the mm emission peak. Interestingly none of these species 
peak exactly toward the mm continuum, but all peak positions 
are shifted a little bit to the south-east (of the order 1 - 2"), 
largely coi nciding with the cm emission peak attributed to the 
thermal jet (Zapata et al. 2005). Although the shift in the spec- 
tral line maps is smaller than the angular resolution of the syn- 
thesized beam, we consider it a real observational feature be- 
cause (a) the line and continuum images are extracted from the 
same dataset, and (b) by fitting peak positions one can achieve a 
significant higher angular resolution than the actual beam size 
HPBW depending on the signal to noise ration S/N (down to 
0.5 X HPBW/(S/N), Reid_et_al, 1988) . A close inspection of 
the higher density tracing '^CO blue and red outflow emission 
shows that both line wings shows strong emission toward the 
south-east of the mm peak as well, and are thus also associ- 
ated with the cm jet. A spatial offset between most molecu- 
lar species - particularly high-density tracers like CH3CN or 
HCOOCH3 - and the mm continuum emission is surprising, 
and various potential origins will be discussed in ^ 



3.3. The molecular outflows 

Figure 13 shows the SMA-only '^CO, '^CO and C"*0 spec- 
tra vector-averaged over all baselines. The '^CO spectrum 
shows a dip around the systemic velocity visr of 59.1 km s"', 
which is not real absorption but caused by the missing short- 
spacing information in the interferometer data. The interfero- 
metric spatial filtering affects the most abundant isotopologue 
'^CO strongest because it is optically thick at the visr, trac- 
ing the extended cloud envelope. In contrast to that, the less 
abundant isotopologues '^CO and C''^0 have lower opacities, 
they are tracing more compact structures and hence do not 
show the pronounced dip at the visr. Nevertheless, already these 
vector-averaged spectra in the uv-domain show '^CO high- 
velocity line wings due to the molecular outflow between 30 
and lOOkms L Although weaker and spanning a narrower ve- 
locity range, the line-wing emission is also found in the '-^CO 
data but barely discernable in the rare isotopologue C"*0. 
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Fig. 4. Integrated images of the observed molecular species. The full and dashed contours present the integrated emission as 
labeled in '^CO and '^CO image. The full and dashed contours in the remaining panels present the positive and negative (due to 
the side-lobe emission and inadequate cleaning) integrated emission as labeled in the panels. The contour levels range from 10 
to 90% (step 10%) of the peak emission reported in Table|3] There we also list the Icr rms of each integrated image. The star 
marks the 1.3 mm continuum peak position from Fig.|2l the open and full triangle mark the cm emission sources with positive 
and negative spectral index, respectively tZapata et al..2005ii) . The synthesized beams are shown at the top-right of each panel. 



Since we have the corresponding CO(2- l) single-dish 
observ ations from the IRAM 30 m telescope jBeuther et al.l 
l2002bl and Fig. |H1 top-left panel), we were able to combine in- 
terferometer and single-dish data and thus recover the miss- 
ing short spacings. The IRAM 30 m data were less sensitive 
than the new SMA observations detecting high-velocity emis- 
sion only approximately between 50 and 70 km s ' . Therefore, 
we confined the SMA data to the range between 47.5 and 



72.5 km s"' to combine both datasets. Figure |6l presents a 
'^CO(2-l) channel map showing the merged SMA+30m data 
between 50 and 70 km s ' and the SMA only data at higher ve- 
locities relative to the visr. Figure|8lshows the resulting outflow 
images covering all spatial scales from the single-dish to the 
interferometer observations. 

The interferometer data alone exhibit CO emission approx- 
imately between 30 and 100 kms"' (Figs.|6l&Q. We tried to 
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Fig. 6. '^CO(2-l) channel map with 5kms ' resolution. The channels between 35 to 45 kms ' and between 75 to 90kms ' 
are produced from SMA only data, whereas the channels inbetween are from the merged SMA+30M dataset (labeled in each 
panel with the corresponding central velocity). For the SMA only data, we plot the contour levels from +0.5 Jybeam 'kms~^ 
in +0.5 Jy beam 'kms ' steps. The Icr rms of the weak emission channel at 90 kms ' is 0.15 Jy beam^'kms The SMA+30m 
data are contoured differently in each panel to account for the dynamic range differences. The start- and step-values are the same 
in units of Jy beam^'kms"'; 50kms"': +1.6; 55kms"': +4.4; 60 kms"': +6.6; 65 kms"': +3.4; 70kms"': +1.0. The crosses 
mark the position of the m m conti nuum p eak, and the open and full triangles mark the positions of the cm sources with positive 
and negative spectral index dZaoata et alJ2005l) . 



Table 2. Observed lines 



Table 3. Line image parameters to Fig.|3 



V 


line 


^lower 


[GHz] 




[K] 


LSB 


219.560 


C"*0(2- 1) 


5.3 


219.798 


HNCO(10o,io - 9o,9) 


48 


219.909 


H^3cO(3i,2-2,,,) 


22 


219.949 


SO(65 - 54) 


24 


220.079 


CH30H(8o.8 - 7i.6) 


85 


220.167 


HCOOCH3(172.i5 - I64.2) 


93 


220.178 


CH2CO(li,„-10uo) 


66 


220.190 


HCOOCH3( 174.13 - I64.12) 


93 


220.399 


"CO{2 - 1) 


5.3 


220.585 


HNCO(10i.9 - 9i.8) 


91 


220.594 


CH3CN(126- lis) 


315 


220.641 


CH3CN(125 - II5) 


237 


220.679 


CH3CN(124- II4) 


173 


220.709 


CH3CN(123 - II3) 


123 


220.730 


CH3CN(122 - II2) 


87 


220.743 


CH3CN(12i - 111) 


65 


220.747 


CH3CN(12o - llo) 


58 


USB 


229.589 


CH30H(154,u - 163,i3)E 


362 


229.590 


HCOOCH3(193.i6- 184,i5)E 


106 


229.759 


CH30H(8i.8-7o.7)E 


77 


229.858 


'"802(42.2 - 3, ,3) 


7.7 


229.864 


CH30H(195.i5 - 204.i6)A+ 


568 


229.939 


CH30H(195,i4 - 204,i7)A- 


568 


230.027 


CH30H(32.2 - 4,,4)E 


28 


230.538 


'2CO(2 - 1) 


5.5 



image the 24 kms"' feature (Fig.|5} but could not identify any 
clear signature associated with IRAS 18182-1433. Therefore, 
we conclude that this spectral feature has to be associated with 
an unrelated fore- or background cloud. The red-shifted feature 
in the SMA-only CO outflow image in Figure |8l (top-middle 
panel) at Declination ~ -11" and R.A. between approximately 
-5" and -20" is an artifact of the image deconvolution and the 
strong side-lobes caused by the insufficient uv-coverage and 
missing short-spacings of the SMA data alone. After adding in 
the short spacing information from the IRAM 30 m observa- 
tions, this feature is largely gone, whereas the southern edge of 
the red-shifted CO cone at Declination 0" is recovered much 
better (Figs.gl&ill. 

We find red and blue outflow emission between 30 and 
70 kms"' south-east of the mm peak, whereas toward the 



line 


Icr 

r mjy km -| 
beam s 


S peak 
r Jy kmi 
beam s 


'^CO(2 - 1) blue 


1040 


37.7 


'2CO(2 - 1) red 


1088 


29.2 


"CO(2 - 1) blue 


192 


4.2 


"CO(2 - 1) red 


224 


5.1 


C"*0(2 - 1) 


440 


13.3 


CH3CN(125 - II5) 


200 


4.6 


SO(65 - 54) 


630 


22.0 


CH30H(8i.8 - 7o.7)E 


480 


24.9 


HNCO(10o.io-9o.9) 


240 


2.4 


H»CO(3,.2-2,.i) 


168 


2.8 


HCOOCH3( 193.16 - 184.i5)E 


224 


4.3 




Velocity (km/s) 



Fig. 5. SMA-only '^CO, '^CO and C'^0(2-l) spectra vector- 
averaged over all baselines. The spectral resolution is the nom- 
inal resolution of 1.1 kms"'. The vis, at 59.1 kms"' is marked 
by a vertical line. 



north-west only red CO emission is detected. The north- 
western emission exhibits a clear cone-like morphology, 
whereas the south-eastern emission appears jet-like with the 
compact component being ass ociated with the the rmal jet fea- 
ture observed in cm emission dZapata et al.lbOOsI) . Integrating 
over slightly different velocity regimes for the blue wing, the 
observed structures do vary a bit: Employing the originally 
from the single-dish data alone identified blue range from 53 
to 57 kms"', the south-eastern jet-feature dominates the blue 
wing of this outflow. However, integrating a little bit broader 
regime from 48 to 58 kms"', one can tentatively identify as- 
sociated structures at offset (10"/ - 13") and (12"/3") which 
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Fig. 7. SMA-only '^CO, '^CO and C'^0(2-l) spectra extracted 
from the data cube after imaging with a spectral resolution of 
5 km s ' . The spectra are averaged spatially over the whole re- 
gion covered by the CO outflow (see top-left panel in Fig.0}. 
The visi- at 59. 1 km s"' is marked by a vertical line. 

Table 4. Line image parameters to Fig.|S] 



Contours 
[Jy beam"' km s"'] 



30m blue 
30m red 
SMA blue 
SMA red 
SMA+30m blue 
SMA+30m red 



66(22)198 
41(21)185 
1.36(1.36)12.24 
2.07(2.07)18.63 
10.92(3.60)32.76 
8.82(441)39.69 



could be part of a blue counter-cone. This potential blue cone- 
like morphology is less prominent than for the red north- 
western wing, and with this dataset we cannot unambiguously 
decide whether the blue outflow wing only shows jet-like emis- 
sion or whether it actually features a cone-like morphology. 

The single-dish data alone indicated an outflow oriented 
roughly in north-south direction whereas the interferometer 
data suggest an orientation approximately in north-west south- 
east direction. The combined data now show that we are likely 
dealing with a system of multiple outflows: The one outflow 
(a) is the north-west south-east outflow revealed in detail by 
the SMA observations. However, most of the single-dish flux 
at blue-shifted velocities is to the south of the mm continuum 
peak, and Figure |S] shows additional red-shifted emission to- 
ward the north-east of it. One possibility would be that the 
whole blue/red outflow structure is due to one inhomogeneous 
wide-angle outflow with its symmetry axis along the jet direc- 
tion. The features observed with the SMA would then poten- 
tially trace the denser parts of the molecular outflow. In con- 
trast to this picture, the cone-like red-shifted emission as well 
as the jet-like blue-shifted features are strongly suggestive of 
at least one more homogeneous outflow in north-west south- 
east direction as discussed above. To account for the additional 
blue-shifted emission in the south and the red-shifted emission 
toward the north-east, we interpret these features as due to a 
second outflow (b) emanating from the same mm continuum 



core. The two outflows are inclined approximately 90° to each 
other (projected on the plane of the sky). The corresponding 
blue and red outflow lobes are for both outflows not exactly 
in opposite directions (as sketched in Fig.|8|i, but appear to be 
bend a little bit. Since we expect that the two outflows are in- 
teracting with each other, such an observed bending is no sur- 
prise but can be attributed to overlapping contributions from all 
outflows. In summary, the merged CO(2-l) dataset indicate a 
quadrupolar outflow system emanating from the main mm con- 
tinuum source in IRAS 18182-1433. 

To increase the complexity of the outflow system from 
IRAS 18182-1433, Fig. |51 presents an overlay of the CO(2- 
1) data from the SMA-i-30m observations with VLA SiO(l- 
0) observations integrated over the velocity range from 40 to 
90kms ' (kindly provided by Yuan Chen et al., in prep.). The 
south-eastern SiO feature is consistent with the blue jet-like 
CO emission, whereas the SiO emission north of the mm con- 
tinuum peak shows a different morphology compared to the 
red CO cone-like outflow. The strong northern SiO peak posi- 
tion appears to be correlated with the northern edge of the CO 
cone in that direction. However, the rest of the SiO emission 
continues even further north, unrelated to the CO outflow. We 
detect additional weaker SiO emission at an offset (-16"/!") 
spatially associated with the western end of the CO cone. The 
additional SiO peak at (-16"/!") and the northern peak corre- 
lated with the northern CO cavity wall allow the interpretation 
that the SiO emission may be due to the same outflow observed 
in CO, just more strongly excited in the cavity walls. However, 
the SiO extension further north is suggestive of an additional 
third outflow so far not detected in the CO emission. 

Figure [TO] presents position-velocity (pv) diagrams of the 
molecular gas along the directions of the two cavity walls of the 
north-western outflow cone. Both diagrams show that the high- 
velocity gas on the blue-shifted south-eastern side of the out- 
flow remains very close to the outflow center whereas the red- 
shifted north-western side resembles more the typical Hubble- 
law with increas ing velocity at fu rther distances from the out- 
flow center (e.g.. lArce et alJ2006l) . This very different blue/red 
outflow pv-characteristics are likely due to the fact that the blue 
CO outflow emission apparently traces predominantly the jet- 
like component also visible in cm and SiO emission, whereas 
the red CO emission is found mainly toward the cavity-like 
walls of the north-western outflow component. 

The combined CO(2-l) data can also be used to esti- 
mate the contributions of the two outflows emanating from 
IRAS 18182-1433. To calculate the outflow mass M (Mbiue, 
Med, Mot) the mom entum p and the energy E, we follow the 
approach outlined in lCabrit & Bertout'('l990') and already em- 
ployed for the single-dish data alone by Beuther et al. (20023)- 
In the following section ( ^3.4> . we will discuss in more de- 
tail the varying CO line opacities, however, for the total pa- 
rameters M/p/E assuming an average '^CO/'-'CO (2-1) line 
ratio of 10 to calculate the opacity-corrected column densi- 
ties of the blue- and red-shifted outflow emi ssion is a reason- 
able approach (Choi et al. 1 993t iLevreaull 1 988: Beuther et alJ 
2002 bl). Furthemiore, \y e use an Ha/'^CO ratio of 89 x 10"^ 
( Cabr it & Bertou3 Il992h . an average temperature within the 
outflow of 30 K, and the velocity ranges of [53,57] and 




Fig. 8. '^CO(2-l) outflow images observed with the IRAM 30 m telescope jBeuther et al."2002b') and the SMA. The top-left 
panel shows the 30 m data only, the top-middle presents the SMA data only. The top-right and the bottom panels then show the 
combined SMAH-30m data. The velocity-integration regimes for the top-row panels and the bottom-right panel are: blue (full 
contours) [53,57] kms"' and red (dashed contours) [62,71]kms"'. The bottom-left panel has a slightly broader blue velocity 
range from [48,581 km s~'. The grey-scale in the top-left panel presents the 1.2 mm continuum emission from the IRAM 30 m 
teeuther et alJ2002al) . The stars in the other panels mark the mm continuum peak from the SMA observations. The open and full 
triangles in the top-middle and bottom panels mark the positions of the cm emission peaks with positive and negative spectral 
index, respectively (Zapata et al. 2005). The (synthesized) beams are shown at the top-left of each panel, and the axes are off'sets 
in arcseconds from the phase center. The hnes in the top-middle panel show the pv-cuts presented in Fig[TO| and the arrows in the 
bottom panel mark the approximate directions of the outflows discussed in the main text. The contour levels are listed in Table|3 



[62,71]kms ' for the integrated blue and red outflow wings. Table 5. Outflow parameter 
respectively. The derived masses, momenta and energies are 
shown in Table |5l 



Although we recover approximately 88% of the red and 
blue wing emission when integrating over exactly the same re- 
gion (see ^2.2> . the derived parameters from the merged dataset 
are a little bit lower because the area used to derive the out- 
flow parameters from the single-dish data alone is larger than 
for the merged dataset (see Beuther et al. 2002b and Fig. |8}. 
Interestingly, both outflows contribute approximately the same 
mass, momentum and energy to the total budget of the whole 
quadrupolar outflow system. 
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Fig. 9. The thin full and dashed contours show the same blue 
and red CO(2-l) observations merged from the SMA and 
IRAM30m data as presented in Fig.|Sl The grey-scale shows 
the VLA SiO(l-O) emission integrated from 40 to 90kms"' 
(kindly provided by Yuan Chen et al., in prep.). The con- 
tour levels of the CO images are presented in Table |3 and 
the SiO data are contoured from 18 to 342 Jy beam" 'km s"' 
(step 27 Jy beam 'km s '). The CO and SiO synthesized beams 
are shown at the bottom-left and bottom-right, respectively. 
The cross marks the position of the mm continuum peak, 
and the open and full triangles mark the positions of the cm 
sources with positiv e and negative spectral index, respectively 
jZapata et al.l2005h . 



3.4. Outflow line opacities 

The simultaneous observations of the '^CO line and their iso- 
topologues '^CO and C'^O (Figs.|5l&0 should in principle 
allow to study the line wing ratios and the corresponding CO 
line win g opacities with c hanging velocity. The '^CO/'^CO 
study bv lChoi et all ( Il993h indicated that the '^CO line opac- 
ities likely decrease with increasing velocity with respect to 
the visr but the data were not sufficient to test this suggestion 
in more detail. Our new data toward IRAS 18182-1433 have 
the advantage that the simultaneous observation of the three 
lines minimizes the calibration uncertainties. However, the dis- 
advantage is that the spatial filtering effect of an interferometer 
can vary for the different isotopologues. This is due to a number 
of reasons, two of them are: (1) Varying line opacities between 
the isotopologues and thus different spatial structures traced 
by each of them. (2) The synthesized maps are affected by 
the side-lobes of the strongest emission features which causes 
different distortion of the weaker emission in various isotopo- 
logues. Since we have no corresponding single-dish '^CO(2-l) 
and C"^0(2-l) observation s, it is hard to quantify t his effect ex- 
actly for all three lines, but lBeuther et alJ have shown 
for SMA observations of Orion-KL that this effect can signif- 
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Fig. 10. Position-velocity diagrams for the '^CO(2-l) SMA 
only outflow observations with a velocity resolution of 
2 km s ' . The top panel is a pv-plot in R.A. at a Declination off- 
set +0.6" (corresponding to the strongest features in the west). 
The bottom panel presents a pv-diagram in north-west south- 
eastern direction with a RA of 38° from north (the cuts are 
marked in Fig|S| top-middle panel). The contour levels are from 
5 to 95% from the peak emission (4.8 Jy beam"' in both pan- 
els). The visr at 59.1 kms"' and the central positions are marked 
by vertical and horizontal lines.. 



icantly change the observed line ratios. Therefore, we stress 
that the absolute line ratios based on the SMA-data alone are 
not reliable. However, even the relative changes in the '2C0(2- 
l)/'''CO(2-l) line ratio is interesting and allows conclusions 
about the changing line-wing opacities. Figure^Jprcsents the 
'^CO(2-l)/'^CO(2-l) line ratio averaged over the whole re- 
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gion of the molecular outflow. We clearly see the trend of in- 
creasing line ratio with increasing velocity. While the ^^CO(2- 
l)/'"'CO(2-l) line ratio is low close to the visr, it increases to 
higher values at higher relative velocities. Although the abso- 
lute numbers suffer from the above discussed spatial filtering 
problems, it appears that the average line wing ratio of 10 sug- 
gested by the observations of .Choi et aL (.1993.) holds for this 
observations as well. 
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Fig. 11. Ratios between the '^CO and '^CO 7 = 2-1 spectra 
shown in Fig. The ratio increases with increasing velocity 
of the molecular outflow. We had to truncate the ratio < 45 
and > 70kms"' because we detect no '-^CO signal above the 
noise at higher velocities. The visr at 59.1 kms"' is marked by 
a vertical line. 



3.5. Temperature from CHi,CN{\2k - lU) 

Since we detected 7 lines of the CH3CN(12i: - lU) A;-ladder 
with k - 0...6, we can utilize the varying excitation levels 
of the Unes with lower level energies Slower between 58 and 
315 K (see table |2j to estimate a temperature for the central 
gas core. We tried to model the whole CH3CN(12j^ - 11^) 
spectrum toward the mm continuum peak position in the lo- 
cal thermodynamic equilibrium using the XCLASS superset to 
the CLASS software developed by Peter Schilke (priv. comm.). 
This soft ware package uses the line catalogs from JPL and 
CDMS (Pov nter & Pick ett 1985': lMuller et amool . Figure[l2l 
shows the observed CH3CN(12t - \ \k) spectrum toward the 
mm peak position and a model spectrum with a temperature 
of 150 K. This model spectrum fits the observation reasonable 
well. The main difference is that the model spectrum is op- 
tically thin whereas the lower line intensity of the observed 
k-3> line indicates moderate optical depth of the CH3CN lines. 
Nevertheless, the observations clearly show that a hot molecu- 
lar core has formed at the center of IRAS 18182-1433. 
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Fig. 12. The full line shows the CH3CN(12t - lU) spectrum 
toward the mm continuum position. The dotted line presents 
a model spectrum with Trot — 150 K and A^chsCn = 3.5 x 
cm 
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3.6. Disk signatures? 

Although these observations were meant to study the molec- 
ular outflow and thus have not extremely high angular reso- 
lution (3" correspond to approximately 13500AU at the as- 
sumed near kinematic distance), we investigated the typical 
high-density hot-core tracing molecules (CH3CN, HCOOCH3, 
HNCO) for potential disk-signatures in their velocity struc- 
tures. We fitted the peak position of each velocity channel 
separately to ach ieve a higher po sitional accuracy down to 
0.5 HPBW/(S/N) (R eid et al.l ll988) and plotted their positions 
in Figure [O] We cannot identify any coherent velocity struc- 
ture in these molecular lines, the velocity distributions appear 
to be random and are off'set from the mm continuum peak. For 
CH3CN and HCOOCH3, the distributions appear somehow in- 
between the main mm continuum peak and the cm feature. 

4. Discussion and Conclusions 

The combined SMA+lRAM30m CO(2-l) and the VLA 
SiO( 1-0) data reveal a multiple outflow system emanating from 
IRAS 18182-1433. One of the two bipolar outflows traced by 
the CO(2-l) emission shows a collimated blue and a cone- 
like red outflow wing. The northern SiO(l-O) emission is sug- 
gestive of an additional third outflow emanating toward the 
north. We do not identify a clear southern bipolar counterpart 
of this potential third outflow. Similar to this, a few other ex- 
amples have recently been reported where the SiO emission 
is tracing only one lobe of an outflow from a massive star- 
forming region (e.g., IRAS 18089-1 732 fBeuflier et alJl2004al 
IRAS 18566+0408, Zhang etal. 2006). Although subject to 
large uncertainties, an analysis of the '^CO(2-l)/'-'CO(2-l) 
line ratios indicates that the CO line wing opacities decrease 
with increasing velocities with respect to the visr. 

All molecular line maps show strong emission approxi- 
mately 1 - 2" south-east of the mm continuum peak, well cor- 
related with a cm emission feat ure maybe due to an underlying 
thermal jet jZapata et al.ll2005h . The offset between the high- 
density tracing molecular lines like CH3CN or HCOOCH3 and 
the 1 .3 mm continuum peak is surprising and mainly two pos- 



Beuther et al.: SMA observations of the HMPO IRAS 18182-1433 



11 




1 0.5 -0 5 -1 1 5 -0.5 -1 1 0.5 -0.5 -1 

Offset RA. ["] Offset. R.A ["] Offset R.A. ["] 



Fig. 13. Positions of different velocity channels for CHsCNCHs - 1 15), HCOOCH3(l 93,16) - 184,15) and HNCO(10o,io - %,9). 
The positions were fitted from the imaged data-cubes with the task IMFIT in MIRIAD. The numbers label the center velocities of 
each channel, and the error-bars are only the statistical errors of the fitting routine. The star marks the mm continuum peak, and 
the 0/0 is the phase reference center. We connected the positions with increasing velocities to stress their random distribution. 



sibilities are likely to explain that observational feature. First, 
the cm emission feature pinpoints the location of a secondary 
source which remains u nresolved in our mm o bservations (see, 
e.g., IRAS 1 629.1-2422. lchandler et 3112001 . This secondary 
source then could be in an evolutionary more evolved state 
altering the chemistry compared to the potentially younger 
source embedded at the location of the mm continuum peak. 
Contrary to that, the cm emission could be solely due to the 
jet itself and not be associated with another proto stellar ob- 
ject (see, e.g., IRAS 16547-4247. lRodriguez et alJ2 00.5). This 
would suggest that the chemistry is strongly altered by the 
molecular outflow, which is particularly surprising for typical 
hot core and dense gas tracers. Only higher angular resolution 
(sub)mm continuum observations will be capable to resolve 
this ambiguity. 

Although our spatial resolution of 3.6" x 2.4" corresponds 
at the assumed near kinematic distance of 4.5 kpc only to 
an average linear resolution of ~ 13500AU, it is interest- 
ing to note that we see no sub-clustering in the 1.3 mm con- 
tinuum data but only one single-peaked s ource. This is sim- 
ilar to the cases of IR AS 20126-H410 4 (Cesar oni et all 1 19991) 
or IRAS 18089-1732 jBeuther et al J 12005b.) . but it is differ- 
ent to the fragmentation of the core in IRAS 19410+2336 
llBeuther & Schilkell2004 . However, one has to keep in mind 
that the observations toward IRAS 19410+2336 had a linear 
spatial resolution of ~ 2000 AU, approximately 6 times better 
than the new data toward IRAS 18182-1433. Therefore, it is 
too early to constrain differences and similarities between the 
fragmentation process es of varying cores, but as pointed out in 
iBeuther et alJ ll200 5b). there are potentially different fragmen- 
tation signatures in young massive gas cores which are worth 
further investigation. The identification of three molecular out- 
flows as well as the two cm emission features observed in the 
vicinity of the mm peak implies the existence of multiple driv- 
ing sources within this central massive star-forming core. It 
would be interesting to identify and characterize the driving 
sources of each outflow independently, but with the given data 
so far we cannot set further constrains on the individual driv- 



ing sources. Again higher angular resolution observations at 
(sub)mm as well as mid-infrared wavelengths are necessary for 
a more accurate characterization of the driving source proper- 
ties. 

Finding the right molecular line transitions to study 
massive accretion disks remains a difficult task. The of- 
ten used molecule CH3CN apparently does not work for 
IRAS 18182-1433, but the recently suggested HCOOCH3 
does not trace any rotating structure in this source as well. 
While in principle it is possible that there is no rotating gas 
structure in IRAS 18182-1433, we consider this as unlikely 
because molecular outflows require accretion and rotating ac- 
cretion disks (e.g., Cesaroni et al. 2006). Therefore, we believe 
that the non-detection of any rotation signature close to the core 
center is mainly an observational problem. Various explana- 
tions are possible in such very young massive and dense molec- 
ular cores, the main four are: (i) The molecular lines are usu- 
ally not only produced in accretion-disk-like rotating structures 
but also in their surrounding cores, which results in various 
overlapping velocity components and no coherent projected 2- 
dimensional velocity signature, (ii) In many cases, the spatial 
resolution is still not sufficient and multiple unresolved sources 
distort the observed velocity field, (iii) Many molecular lines, 
even the high-density tracer, are strongly affected by the molec- 
ular outflow distorting the observed velocity field, (iv) We are 
dealing with extremely high gas column densities and thus po- 
tentially high opacities. This way, one would not even trace 
the inner rotating structure. These observations confirm that the 
quest for the best disk-tracing molecular fine transitions is still 
open. 
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